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Abstract 

We have simulated four-dimensional interacting self-avoiding trails (ISAT) on the hyper- 
cubic lattice with standard interactions at a wide range of temperatures up to length 4096 and 
at some temperatures up to length 16384. The results confirm the earlier prediction (using 
data from a non-standard model at a single temperature) of a collapse phase transition 
occurring at finite temperature. Moreover they are in accord with the phenomenological 
theory originally proposed by Lifshitz, Grosberg and Khokhlov in three dimensions and 
recently given new impetus by its use in the description of simulational results for four- 
dimensional interacting self-avoiding walks (ISAW). In fact, we argue that the available data 
is consistent with the conclusion that the collapse transitions of ISAT and ISAW lie in the 
same universality class, in contradiction with long-standing predictions. We deduce that 
there exists a pseudo-first order transition for ISAT in four dimensions at finite lengths while 
the thermodynamic limit is described by the standard polymer mean-field theory (giving a 
second-order transition) , in contradiction to the prediction that the upper critical dimension 
for ISAT is du = 4. 
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1 Introduction 



The geometric collapse of polymers in dilute solution, being one of the fundamental and most 
well studied phase transitions in statistical mechanics attracts continuing interest from both 
the theoretician attempting to understand the subtleties inherent in this phenomenon and the 
experimentalist wanting a firm foundation for interpreting their results on more complicated 
systems. One approach has been the use of lattice models of various types of self-avoiding paths 
to describe the geometry of linear polymers with added local potentials to broadly account 
for the complex solvent-polymer and polymer-polymer interactions. The canonical lattice path 
utilised in this regard has been the self-avoiding walk (SAW) as it explicitly demonstrates the 
"excluded volume" expected of physical polymers. However, much of our understanding of 
three-dimensional polymer collapse ^, ^ has come from calculations using the continuum 
Edwards model ||^, |5| which is based upon intersecting random paths, rather than interacting 
self-avoiding walks (IS AW). Hence it is of interest to study lattice models closer in nature to 
the Edwards model. One such model is based on the lattice paths known as self-avoiding 
trails (SAT) ^, ^, ^ rather than self-avoiding walks. The type of interaction considered 
(contact versus nearest-neighbour for ISAW) is also more closely analogous to the Edwards 
model. It is also worth noting that the upper critical dimension for polymer collapse is expected 
to be du = 3, and therefore many subtle logarithmic corrections are predicted to occur in 
three dimensions. As is often the case in the study of critical phenomena our understanding is 
enhanced by considering the dimension of the problem as a parameter that can be varied. Any 
discrepancies between the behaviour of competing models in four dimensions for example would 
impinge on our interpretation of results in three dimensions. 

Self-avoiding trails (SAT) or trails for short are paths on a lattice which have no two steps 
on the same bond of that lattice but may occupy the same site. This restriction is sometimes 
referred to as bond-avoiding, in contrast to self-avoiding walks (SAW) which are site-avoiding 
(that is, no two vertices of the walk may occupy the same site on the lattice). Clearly walks 
are, by default, also bond-avoiding. Trails possess an excluded volume effect and it is fairly well 
established that SAT and SAW are in the same universality class |ll| which describes good 
solvent polymers. (For a more recent discussion of the subtle differences between walks and trails 
regarding corrections to scaling see [^.) It has been shown ^ that there should exist a collapse 
transition when contact attraction is added to the trail model: this model of polymer collapse 
is known as interacting self-avoiding trails (ISAT). Moreover, Shapir and Oono have argued 
that this point should be tricritical in nature, as it is at the ISAW collapse point. However, they 
predict that ISAW and ISAT are in different universality classes. Importantly, while the upper 
critical dimension for ISAW is expected to be d„ = 3, the Shapir-Oono field theory gives = 4 
for ISAT. Therefore, this implies generically that logarithmic corrections occur at the 0-point in 
four dimensions, which presumably should occur at a finite temperature. On the other hand the 
above discussion about the hypothesised equivalence of the critical phenomena of self-avoiding 
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walks and the Edwards model may lead one to the opposite conclusion, namely that the collapse 
transitions of ISAT, ISAW and the (3-parameter) Edwards model lie in the same universality 
class. Computer simulations in two and three dimensions have given conflicting results and so 
it is of interest to ascertain whether computer simulations in four dimensions may shed light on 
the question of whether ISAW and ISAT share the same collapse universality class. 

Recent developments have been two-fold. Firstly, preliminary evidence of collapse in any 
polymer model in four dimensions was presented in using the so-called kinetic growth trail 
which effectively simulates ISAT with a particular fixed set of non-standard Boltzmann weights 
(and so fixed temperature). It was argued from the results of the simulations that this set 
of Boltzmann weights corresponded to a so-called ^-point where the ISAT behaves in a way 
predicted, by the mean-field theory of the polymer collapse, to be precisely at the collapse 
transition point. Plausibility arguments then deduced the existence of the collapse transition 
as the temperature is varied through this value. Secondly, and more startling have been the 
results of simulations [ll4| , ^] of interacting self-avoiding walks in four dimensions. A dramatic 
collapse transition was located in this model using simulations up to length 32678. The nature 
of this transition was described as a pseudo-first order transition because at any finite length the 
distribution of the internal energy was bimodal when the temperature was near that at which 
the specific heat attained its maximum (one indicator of the transition's position), while it was 
argued that the effective latent heat would slowly decay to zero in the thermodynamic limit (as 
the length diverges) with an anomalous exponent. On the other hand the specific heat seemed 
to diverge as the length was increased. Despite these peculiar findings the simulational results 
were able to be interpreted in terms of a framework explained some time ago by Khokhlov [p!^ ], 
who applied the non-standard mean-field approach of Lifshitz, Grosberg and Khokhlov (LGK) 



1 17, |l^, 19 1 to arbitrary dimensions. This theory is based on a phenomenological free energy 
in which the competition between a bulk free energy of a dense globule and its surface tension 
drive the transition. The consequences of this surface free energy were largely ignored in the 
polymer literature until recently, when its effect on the scaling form of the finite-size partition 
function was proposed and confirmed |20, 21, 22, 23]. While markedly different for finite lengths 
this theory still predicts a standard Gaussian 0-point and in the infinite length limit a second 
order phase transition with a jump in the specific heat (no divergence as the temperature is 
varied: the divergence as length is varied exists but "disappears" in the thermodynamic limit). 

In this work we have simulated the standard ISAT model on the four-dimensional hypercubic 
lattice over a range of temperatures using the PERM algorithm [p4| in a similar fashion to sim- 



ulations implemented for ISAW |14]. This algorithm is particularly efficient in high dimensions 
and so is well suited to these simulations. We find evidence that parallels the findings in the 
ISAW simulations |14] for a collapse transition with the signature of pseudo- first order type. 
The transition is, if anything more pronounced, with an equivalent strength at shorter lengths. 
We demonstrate that our data is at least as, if not more, consistent with the LGK theory as 



is the ISAW results [14|. Hence we deduce that the collapse transitions of ISAW and ISAT in 
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four dimensions are in the same universality class. It follows from these conclusions that the 
upper critical dimension for collapse in ISAT is du < 4, and that possibly = 3 in accord 
with the standard polymer theory and with the logarithmically modified Gaussian state found 
in three-dimensional kinetic growth trail simulations |25|. 

2 The ISAT model and scaling theories 
2.1 The ISAT model 

We define the ISAT model on the four-dimensional hypercubic lattice in the following way. The 
lattice has coordination number 8 and we consider configurations of trails, or bond-avoiding 
walks, of length (bonds) starting from a fixed origin. Let rrik, k = 1, . . . ,4 be the number 
of sites of the lattice that has been visited k times by the trail so that ^ km^ = A + 1. The 
partition function of a very general interacting trail model is 



ZNiuJ2,u;3,u;^) =Y,^r^r<' , (2-1) 



'Pn 

where oj^ is the Boltzmann weight associated with A;- visited sites. The canonical model is one 
where every segment of the trail at some contact site interacts with every other segment at that 
site, so that 

(jfc = u;(2) for /c = 2,3,4 , (2.2) 
with Lo = i02- This implies that in our specific case 

UJ2 = UJ, U!s=UJ^, UJ4=UI^. (2-3) 

The Boltzmann weight co = e^^ is associated with a energy of contact — e so that to > 1.0 for 
attractive interactions. As we are only interested in the attractive case here one can set e = 1 
for convenience. We define a reduced finite-size free energy per step K]\[{uj) as 

kn{oj) = ^ log Zn{uj). (2.4) 

The usual free energy is related to this by —(3Fn = Akjv(<^). 

The average of any quantity Q over the ensemble set of allowed paths 17 at of length A is 
given generically by 

{Q)n{u:) = (2.5) 
where m = m2 + 3m,3 + 6777,4. We define a normalised finite-size internal energy per step by 

U^iio) = ^ , (2.6) 



and a normalised finite-size specific heat per step by 

C.(.) = <!=!>^ . (2.7) 

These quantities are related in the usual way to the reduced free energy via C/jv = dK^/dloguj 
and Cat = dU^ /dlogu. Note that because of the absent factors of temperature both Un and 
Civ can take on non-zero values for u = 1.0. 

The thermodynamic limit in this problem is given by the limit iV — > oo so that the thermo- 
dynamic free energy per step foo{i^) is given by 

-/3/oo('^) = Hooii^) = lim KAr(u;) . (2.8) 

AT— >oo 

This quantity determines the partition function asymptotics, i.e. Zn{u!) grows to leading order 
exponentially as with /x(a;) = e'*°°^'^). 

In our simulations we calculated two measures of the polymer's average size. Firstly, specify- 
ing a trail by the sequence of position vectors tq, ri, r^r the average mean-square end-to-end 
distance is 

{RDn = {{tn - ro) • {VN - ro)) • (2.9) 
We shall use the symbol jy to be equivalent to 

Rl^ico) ^ {RDn. (2.10) 
The mean-square distance of a site occupied by the trail to the endpoint, fq, is given by 

1 ^ 

{RDn = i;((r. - ro) ■ (r, - ro)) . (2.11) 

1=0 

Again we define 



<iv(^) ^ {Rin)N. (2.12) 

We also define the ratio 

Bn{u) = (2.13) 
which should have a universal limit in each critical phase of the model. 



2.2 Standard polymer scaling theory 

Let us expand our view in this section to general polymer scaling theory. We first assume that 
there is a single collapse transition at some value of temperature and explore the four-dimensional 
behaviour we might expect from the above defined quantities in each of the phases. The basic 
physics of the coil-globule (collapse) transition can be understood by the consideration of the 
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average size of the polymer, R]\f, either Re,N or Rm,N-, as a function of length in each of the 
phases, so let us consider this first. Let us define an effective Boltzmann weight uj = — of 
course, in any particular model or physical system the associated potentials may be different. 
Generally one always expects that 

i?^ ~ a(w) iV^z. ^ N^oo (2.14) 
for any fixed value of w. In four dimensions at infinite temperature, u; = 1, it has been predicted 



im that 

Rl ^ a+ N {\og{N)Y'\ (2.15) 

If there does exist a collapse transition then one would expect that this scaling extends (with a 
constant a"*" that depends on temperature) down to the transition point. In the collapsed phase 
the polymer is expected to assume a dense configuration on average and hence the globular value 
of the radius-of-gyration exponent is = l/d = 1/4 [|| with 

R% a- {uj) N^/^ . (2.16) 

Finally at some finite transition temperature 1.0 < oot < oo a Gaussian scaling of the radius of 
gyration should occur, that is 

R% ~ N, (2.17) 

so that ut = 1/2. This Gaussian scaling is often used (theoretically at least) to define the 0-point 
UJ = iOg of an isolated polymer so that cot = ujq. The universal ratio i?Ar is expected to converge 
to the value B^o = 1/2 both in the swollen phase and at uj0. However, one would expect slow 
logarithmic corrections for uj < ujg and algebraic corrections at ujq. For uj > ujg the phase is 
no longer expected to be critical and so B^o is no longer universal and may be a non-constant 
function of uj. 

One can also consider the scaling of the partition function in each of the regimes, given that 
there is a transition. For high temperatures 1.0 < u; < ujg one expects the infinite temperature 



behaviour, which is [26| 



Zn ~ b+{u;) /i(cj)^ (logiV)^/^ , (2.18) 



while at low temperatures |20| one expects asymptotics of the form 



Zn ~ b-{uj) fi{uj)^ /is(w)^'^* N3 (2.19) 

where fj,s is related to the surface free energy of the polymer globule and the exponent g need 
not be universal (we only write it for completeness of the asymptotic form). For uj = ujg one 
expects 

Zn ~ K^e)^ (2.20) 



as a reflection of Gaussian behaviour. 

In the thermodynamic hmit the thermodynamic functions fooi^^), Uoo{uj) and Cqo{uj) are all 
expected to be analytic functions of a; except at loq. By using the correspondence to the tricritical 
model m the mean field theory would imply that the specific heat had a jump discontinuity at 
ujg (the associated exponent a = 0). Of course, for finite N there is no sharp transition for an 
isolated polymer (unless one examines a macroscopic number of such polymers). 



2.3 LGK theory 

We now provide a brief review of the predictions of the theory of Lifshitz, Grosberg and Khokhlov 
(LGK) 1 19 1 as applied to four-dimensional polymer collapse by Khokhlov |1(:]. Firstly, there 
exists a state where the excluded volume property of long chain molecules is exactly cancelled 
by the attractive interactions between parts of the polymer as mediated by the solvent. This is 
the 0-state. Secondly, when the attraction becomes even stronger there eventuates a globular 
state where the polymer behaves as a liquid drop. The results of the theory are based on 
a phenomenological free energy of that globular state relative to the free energy of the pure 
Gaussian state of the 0-point at Tg. Hence the condition applied to find the finite-size position 
of the transition is to equate the relative free energy to zero. The relative free energy is given 
as a sum of bulk and surface contributions which are, in turn, given in terms of the second and 
third virial coefficients, the length of the chains, and the linear size of the polymer found from 
the globular density. In particular both the bulk and surface free energies are proportional to 
the square of the second virial coefficient. It is assumed that on approaching the 0-point the 
second virial coefficient goes to zero linearly with temperature while the third virial coefficient 
remains non-zero. Note that this implies a quadratic dependence of the bulk free energy on 
the distance to the 0-point. Since the free energy has exponent 2 — a this implies an exponent 
a = (assuming that this part of the free energy is singular). Therefore a second-order phase 
transition occurs in the thermodynamic limit. 

It is further assumed that the density in the globule is proportional to the second virial 
coefficient and hence also goes to zero linearly with temperature on approaching the 0-point 
{P = 1). Again using an effective Boltzmann weight u = and defining the transition as when 
the free energy is zero, Khokhlov ||l^ finds a finite-size transition temperature^ uJc,n which 
approaches the ^-temperature (cog) as 

s 

for some constant s. That is, the polymer collapse shift exponent is 1/3. The width of the 

transition region Auj at finite N can be found from the free energy rewritten in terms of this 

^for the sake of ease of expression in this section we will use the word "temperature" to mean the effective 
Boltzmann weight 



transition temperature to scale as 



Ao; ~ ^ (2.22) 



for some constant w. That is, the polymer collapse crossover exponent is 2/3. Hence note that 
the size of the crossover region is asymptotically small relative to the shift of the transition. 



Following the work of Lifshitz, Grosberg and Khokhlov |1S] one can also calculate the 
change in the internal energy over the crossover width of the transition Auj as the latent heat 
(or "heat of the transition") by using expression of the free energy in terms of the transition 
temperature as 

^^-J^s- (2-23) 
The corresponding height of the peak in the specific heat is 

CTvKiv) ~ /i" iV'/^. (2.24) 

So to summarise the LGK picture, the theory predicts a thermodynamic second-order tran- 
sition at a Gaussian ^-point with a jump in the specific heat. For finite polymer length this 
transition is shifted below the 0-point by a temperature of the order of 0{N^^/^) with the 
width of the transition of the order of 0{N^'^^^). Over this width there is a rapid change in 
the internal energy that scales as 0{N^^^^): the important point here of course is that this 
tends to zero for infinite length so the effect of the peak in the specific heat is scaled away for 
N large, leaving a finite jump in the thermodynamic limit. To understand this further let us 
consider the distribution of internal energy as a function of temperature and length. For any uj 
below ujg and well above uJc,n one expects the distribution of internal energy to look like a single 
peaked distribution centred close to the thermodynamic limit value: a Gaussian distribution 
is expected around the peak with variance 0(A^~^/^). In fact, this picture should be valid for 
all temperatures outside the range [uJc^n — 0{N^'^/'^),uJc^i^ + 0{N^'^/'^)\. When we enter this 
region we expect to see a double peaked distribution as in a first-order transition region. For 
any temperature in this region there should be two peaks in the internal energy distribution 
separated by a gap 5U of the order of 5U ^ AU cx 0{N~^/^). Each peak is of Gaussian type 
with individual variances again of the order of 0{N^^^'^). Hence as N increases the peaks will 
become more and more distinct and relatively sharper but the peak positions will be getting 
closer together. We refer to this scenario as a pseudo-first-order transition or, more correctly, 
as first-order-like finite-size corrections to a second-order phase transition. If there were a real 
first-order transition then the distance between the peaks should converge to a non-zero con- 
stant. On the other hand the transition is not a conventional second-order phase transition with 
a well defined limit distribution of the internal energy that is simply bimodal. 
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3 Simulational results and analysis 



We have simulated ISAT on a 4-dimensional hyper-cubic lattice using the Pruned-Enriched 
Rosenbluth Method (PERM), a clever generalisation of a simple kinetic growth algorithm [ p^ , 
p7| . PERM builds upon the Rosenbluth- Rosenbluth method [28|, in which trails are generated 
by simply growing an existing trail kinetically but overcomes the exponential "attrition" and 
re- weighting needed in this approach by a combination of enrichment and pruning strategies. 
Our implementation here follows our previous IS AW work Briefly, we chose upper 

and lower thresholds and W'', for enrichment and pruning respectively, proportional to the 
current estimate of the average weight of a trail at length N, {Z]\f)/sjsr, where sat is the number 
of generated samples at length N, and {Zi\f) is the current estimate of the partition function at 
length N. That is to say, VF^ = c^{Z]y) /s^, W^f = c^^q {Z jq) / s n ■ In order to enforce an even 
sample size distribution we allowed for dynamic adjustment of and c^, keeping the quotient 
of the thresholds Q = c^/c^ constant. To stabilise the dynamic adjustment, we enforced 

> c^in ^ ^^max- ™ the ISAW work, we chose c^^^ = 2 and c^rnax — 1/2- For each 

run, we attempted to choose the smallest threshold quotient Q for which we could obtain an 
even sample size distribution. 

Each run had a maximum length N^ax set and while individual runs gave information about 
shorter lengths we collected data from independent runs at some shorter lengths to guarantee 
statistical independence. Simulations were conducted with the maximum lengths N^ax set 
to 512, 1024, 2048, and 4096, with values of uj ranging from 1.0 to 2.07 for Nmax = 512, 
from 1.0 to 2.00 for Nmax = 1024, from 1.0 to 1.78 for Nmax = 2048, from 1.0 to 1.67 for 
^max = 4096. We also ran many closer spaced simulations in the range of uj from 1.4 to 1.42 
at length Nmax = 16384. At each fixed oj, we generated at least 10^ trails. To illustrate the 
computational effort, the generation of a sample of size lO"^ at length Nmax = 16384 took about 
2 weeks CPU time on a 600 MHz DEC Alpha. The threshold quotient Q used ranged from 10 
to 80 with larger values of Q needed for higher lo. We also performed one large-scale simulation 
deep in the collapsed regime, with Nmax = 512 and uj = 4.0, for which it was necessary to 
increase the threshold quotient up to Q = 1000. 

We computed statistics for -Rg^ and RmN' partition function Zjy, the internal energy 
Un and specific heat Cat. Moreover, we generated the distribution of the number of interactions 
at Nmax- The distributions obtained at various temperatures were then combined using the 



multiple histogram method |2£]. 

The disadvantage of PERM is that due to the enrichment the generated data is not inde- 
pendent. All the data generated during one "tour", i.e. between two successive returns of the 
algorithm to length 0, is correlated. Therefore, we kept track of the statistics of tour sizes t 
to get a rough idea of the quality of the data. In our statistical evaluation we use (somewhat 
arbitrarily) the quotient of sn and \/Jt^ as a measure of an effective independent sample size. 
This is correct as long as the tour sizes don't fluctuate too strongly, and, more importantly. 
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as long as individual tours explore the sample space evenly. When simulating in the collapsed 
phase, both of these assumptions break down, and the sample is dominated by few huge tours. 
Moreover, the pruning and enrichment rates become so large that the efficiency of the algorithm 
is significantly decreased. Error bars are only given for high temperature and ^-point simulation 
figures (figures |^ and |^ and are based on the method described above, although we always 
computed error estimates. No error bars are given in the rest of the figures because of the 
performance of the algorithm in and near the collapsed phase, even though the data seemed 
converged sufficiently. 

Let us first discuss the scaling of the mean-squared end-to-end distance normalised by trail 
length, iig n/^ ^'^'i mean-squared distance of a site occupied by the trail to the trail's end- 
point, Ft^ n/^ ■ swollen phase, our results are in correspondence with the logarithmic 
corrections seen by Grassberger As in that paper, we observe that -Rg^ grows faster than 
N for oj near 1.0, and fitting to A^(logA^)'^ at = 1.0 gives an effective exponent close to that 
predicted by field theory (1/4). This value shifts as u) is increased indicating the presence of 
strong temperature-dependent correction terms. 

By considering when the quantity n/^ approaches a constant we narrowed our search 
for the ^-point to the region uj = 1.40 to w = 1.42. In this region we extended our simulations 
to trails of length 16384. Figure |l] shows a plot oi j^/N versus 1/A^ for values of uj between 
1.408 and 1.420. At uj = 1.414(3) we have an approximate linear asymptotic dependence of R^ ^ 
on N . Moreover, at a; = 1.414 we estimate from our data Bn = R^ n/^c N ~ 0.5000(2), which 
is also indicative of Gaussian behaviour: the precision of this estimate stems from the weakness 
of the corrections to scaling at this point. 

As shown in Figure ^ for uj = 4.0, R^ jy changes non-monotonically in A^! After an initial 
increase, the size of the polymer actually starts to shrink around A^ = 50 as it undergoes collapse 
corresponding to a rapid increase of the density. For large enough A^, we expect to see the true 
collapsed behaviour, i.e. R^ growing again as N^^'^, but while we see R^ jy just starting to 
increase again, the asymptotic regime is beyond the reach of our PERM simulations on current 
computer hardware. 

Let us now discuss the scaling of the partition function. The swollen phase and the 0-point 
behaviour can also be clearly identified from the free-energy scaling. In the swollen phase we 
find again the same behaviour as |^^. The presence of logarithmic corrections is consistent with 
our data. At a; = 1, we estimate = ^sAT = 6.926080(2). In the ^-region, an analysis shows 
that here scales as jj^ with weak 1/A^ corrections. Figure ^ shows Zjy/Z"^^^ plotted versus 
1/A from which we estimate the 0-point to be uoo = 1.414(3) and fig = 7.0016(6). (At fixed uo, 
the accuracy is of course higher: for uj = 1.414, we estimate fi = 7.0015714(5).) In the collapsed 
region, one expects the finite-size free energy to have a strong correction term of the order N~^^^ 
due to surface effects. Figure |^ shows this for u> = 4.0. As argued above, the globule starts to 
collapse when the length is above A = 50, and we notice here the onset of a corresponding strong 
change in the behaviour of the finite-size free energy around this length (A~^/^ ~ 0.35). Even 



Q 



though we cannot simulate long enough chain lengths to clearly determine the precise nature of 
the correction term, our data is certainly compatible with a N~^/'^ correction for N~^/'^ < 0.3 
(i.e. iV > 150). 

In order to study the collapse transition more closely, we now focus our attention on the 
internal energy and specific heat. As can be seen from Figure |5|, the specific heat has a sharply 
peaked graph for each length that becomes more sharply peaked as increases. The transition 
region becomes sharper and stays well separated from the ^-point, even though the location of 
the transition (peak in the specific heat) approaches the ^-point slowly. The scaling of the shift 
of the transition towards the 0-point, Wc.Af — wg, and the sharpening of the transition width, 
Au!, are both shown in Figure ^. Here, we defined the location of the collapse transition by 
the location of the specific heat peak, and the width of the transition is given by the interval 
in which the specific heat is greater or equal to half the value of the peak height. Expecting 
from the LGK theory that uJc,n — ^e scales as N~^^^ and that Au; scales as N~'^/^ , we plot both 
N'^^'^iiOc^N — ^e) and /S.ujN'^/^ versus N^"^^'^ which was chosen empirically. Both quantities can 
be seen to be asymptotic to constants: on the graph extrapolations give non-zero intercepts. 
Hence, Figure ^ shows that the LGK predictions are compatible with our simulations. We do 
note that the corrections to scaling for Aw are much larger than for uJc,n — ^e- 

The character of the transition becomes apparent if one plots the internal energy density 
distribution (rescaled density of interactions) at the finite-size collapse transition, lOc,n- Figure 
^ shows the emergence of a bimodal distribution. At length 512 one sees a slight non-convexity, 
which at length 4096 has evolved into a distribution dominated by two sharp and well-separated 
peaks. The values of the minima and maxima of the distribution are different by two orders 
of magnitude. This bimodal distribution means that as lo is increased through the transition 
region the density distribution switches from the peak located at a small value of contacts to 
the peak located at a larger value of contacts, corresponding to a sudden change in the internal 
energy. In the collapsed phase, the width of the peak is much wider than in the swollen phase, 
implying a larger specific heat. It is this difference between the swollen and collapsed phases' 
specific heats that will eventually become the thermodynamic second order jump. The rapid 
first-order like switch between two peaks in the distribution becomes more pronounced at larger 
polymer lengths since the depth of the "valley" between the two peaks becomes relatively larger. 

Continuing with the scaling predictions from Khoklov theory, a suitably defined finite-size 
latent heat, AQ, should tend to zero as A^^/^ in the thermodynamic limit. One possible measure 
of this latent heat is given by the product of specific heat peak Cn{'^c,n) and specific heat width 
Au, and another is given by the distance 5U of the peaks in the bimodal internal energy 
distribution. Figure ^ shows the behaviour of both of these quantities. One notices two things 
from this figure. Firstly, it indicates that C]^{uJc^isi)A.uj decreases to zero linearly in N^'^/^ as 
predicted. However, even at length A^ = 2048 (A~^/'^ ~ 0.08) there is considerable discrepancy 
between the two quantities plotted, so that one needs to be cautious in the interpretation of the 
scaling behaviour. The explanation for the discrepancy between the two quantities is of course 
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that in order to observe the asymptotic behaviour the two peaks in the histogram have to be 
well separated and distinct, and that this is only really the case when N is of the order of 10'^. 
We caution that Figure ^ alone is not sufficient to discriminate between the scenario proposed 
here and a real first-order transition in the thermodynamic limit, but we believe the rest of our 
data and other theoretical facts provide a more consistent picture. 



When comparing our data with the simulations for IS AW |15] we note further that the 
bimodal distribution emerges for ISAT at much shorter configurations, so that the peaks in the 
distribution for ISAT at length N = 512 are already more pronounced than the peaks in the 
distribution for ISAW at length N = 2048. To quantify this observation, we turn to the scaling 
predictions of LGK theory. An important parameter in the theory is the quotient a'^/v, where 
a is the mean-square distance between two subsequent monomers (repeated unit element of the 
polymer: equivalent to occupied sites of the lattice model) along a chain and v is the effective 
excluded volume of a monomer, defined via the vanishing of the second virial coefficient at the 



^-temperature. For instance, the shift of the transition temperature (cf. equation 2.21 ) is given 
more explicitly by 

/ ~ 4\ 1/3 

where s is a constant proportional to the quotient of the third virial coeffficient and the excluded 
volume squared. From Figure ^ we estimate that N^^^{uJc^N — '^e) asymptotes to 3.4(1) for ISAT, 
and for ISAW we estimate for the same quantity the value 0.92(3) [^]. Identifying a with the 
lattice constant, which in both models is set equal to one, we can get a rough estimate for the 
relative size of the effective excluded volume v in both models. We obtain 

^ ^ 0.034^ (3.26) 

VSAW SSAW 

and thereby quantify the intuitive notion that the excluded volume effect is numerically "weaker" 
in trails than in walks, though of the same basic type. 

In conclusion, our ISAT simulations elucidate further the structure of the polymer collapse 
transition in four dimensions. We show conclusively that there is indeed a collapse transition 
at a finite temperature. Secondly, we find evidence for a ^-temperature at which the polymer 
is well approximated by Gaussian behaviour as well as for a collapse transition which is well 
separated from the 0-point. The collapse transition shows many first-order like features, such 
as a bimodal distribution in the internal energy. An analysis of the scaling behaviour of this 



transition in the context of the theory of Lifshitz, Grosberg and Khokhlov |19, 16| shows that 
a consistent interpretation of these findings is that of first-order like finite-size corrections to a 
thermodynamic second-order transition. These finding are essentially the same as those made 
recently for ISAW collapse in four dimensions |lj]. Consequently, we deduce that the upper 
critical dimension for ISAT is (i^ < 4 (most likely 3) and not = 4 as was previously predicted. 
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Figure 1: RI^n/N versus 1/N in the 9-region: u = 1.408, 1.411, 1.414, 1.417, 1.420 from top to 
bottom. 



Figure 2: ^ versus N for u = 4.0 up to length 512. 



Figure 3: Zn/Z"^^^ versus 1/N in the 9-region: uj = 1.408, 1.411, 1.414, 1.417, 1.420 from 
bottom to top. 



Figure 4: Finite-size free energy kn versus N for u> = 4.0. 



Figure 5: Specific heat Cn versus uj for lengths 512, 1024, 2048, and 4096 from right to left 
respectively, using the multi-histogram method. 



Figure 6: Scaling of the transition: shift and width of the collapse region. Shown are the scaling 
combinations N^/^{oJc,n — <^e) o.f'd iV^/^Ao; versus N~'^/^. 



Figure 7: Internal energy density distributions at uic,N for 512 and 4096. The more highly peaked 
distribution is associated with length 4096. 



Figure 8: Scaling of the latent heat AC/.- our two measures of AC/, Cn{(^c,n)^^ o-nd peak 
distance 5U are plotted versus N^^/^. 



1/1 



1.36 



RIn/N 1.35 - 



1.34 



i 1 



I 



I 



I 



I 



i 

i 



0.0000 0.0002 0.0004 0.0006 

l/N 



0.0008 



i 



0.0010 



Figure 1: Rl^pf/N versus l/N in the 6-region: lo = 1.408, 1.411, 1.414, 1.417, 1.420 from top to 
bottom. 
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Figure 2: versus N for a; = 4.0 up to length 512. 
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Figure 3: Z^jZ]^!^ versus 1/N in the 9-region: lo = 1.408, 1.411, 1.414, 1.417, 1.420 from 



bottom to top. 
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Figure 4: Finite-size free energy kn versus N for u = 4.0. 
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Figure 5: Specific heat Cn versus to for lengths 512, 1024, 2048, and 4096 from right to left 
respectively, using the multi-histogram method. 
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Figure 6: Scaling of the transition: shift and width of the collapse region. Shown are the scaling 
combinations N^/^{ijJc,n — ^e) o,nd A^^/^Ao; versus N~'^/^. 
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Figure 7: Internal energy density distributions at lUc,n for 512 and 4096. The more highly peaked 
distribution is associated with length 4096. 
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